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1.  Introduction 


In  many  optical  systems  employing  lasers,  it  is  important  to  obtain  the  maximum  far-held  (or 
focal-plane)  irradiance,  for  a  given  output  laser  power,  on  or  near  the  optic  axis.  Examples  of 
such  systems  include  optical  radars,  where  one  wants  to  obtain  the  highest  possible  signal  return, 
and  optical  communications  systems,  where  one  wishes  to  operate  under  conditions  of  the  highest 
possible  signal-to-noise  ratio.  In  these  systems  and  others,  one  is  often  concerned  with  the  far- 
held  radiation  pattern  that  is  well  within  the  central  lobe  of  the  transmitted  beam,  the  correspond¬ 
ing  details  of  the  much  weaker  side  lobes  being  of  less  relative  importance.  Since  all  systems,  in 
practice,  operate  in  the  presence  of  some  (unintentional)  random  angular  jitter,  it  is  of  concern  to 
know  what  the  optimum  truncation  of  the  laser  output  beam  should  be  in  order  to  obtain  the 
highest  possible  average  on-axis,  far-held  irradiance.  Here,  random  angular  jitter  (e.g.,  due  to 
random  vibrations  of  the  optical  platform  mount  and  propagation  through  the  atmosphere  where 
clear  air  tuibulence-induced  wave-front  tilt  produces  eHects  fonnally  equivalent  to  random  plat¬ 
form  jitter^)  refers  to  the  uncertainty  of  the  instantaneous  direction  of  the  optical  axis  with  respect 
to  a  nominal  z-axis  that  is  parallel  to  the  direction  of  the  desired  line  of  sight 

Specihcally,  here  we  model  random  angular  jitter  by  the  commonly  used  assumption  that  it  is 
normally  distributed  about  both  the  x-  and  y-axes  with  zero  mean  and  equal  variances.  As  we 
show  below,  the  far-held  beam  pattern  of  a  truncated  Gaussian  is,  over  the  parameter  range  of 
interest  here,  very  neatly  of  Gaussian  shape.  This,  together  with  the  assumption  of  nonnally  dis¬ 
tributed  jitter  ,  allows  one  to  obtain  an  analytical  expression  for  the  corresponding  average  beam 
irradiance  distribution  in  the  presence  of  jitter.  Maximizing  the  on-axis  mean  irradiance  as  a 
function  of  the  ratio  of  the  beam  diameter  in  the  aperture  to  the  truncation  diameter  yields  the 
optimum  value  of  the  beam  diameter  for  a  given  aperture  diameter. 

In  Section  2,  we  present  a  general  expression  for  the  far-held  irradiance  distribution  in  the 
paraxial  approximation  of  a  circular  iqrertute  illuminated  by  a  Gaussian  beam  in  the  absence  cf 
jitter.  We  then  briehy  review  the  well-known  optimum  truncation  of  the  beam  through  a  hxed 
circular  aperture.  Next,  we  derive  an  accurate  analytical  approximation  within  the  central  lobe  of 
the  exact  irradiance  distribution  over  a  parameter  range  of  interest  Based  on  this  approximation, 
in  Section  3,  we  derive  a  corresponding  analytical  approximation  for  the  mean  irradiance  distri¬ 
bution  in  the  presence  of  random  jitter.  We  then  maximize  the  mean  on-axis  irradiance  as  a 
function  of  the  ratio  of  the  beam  diameter  in  the  aperture  to  the  aperture  diameter  to  obtain  the 
optimum  value  of  the  beam  diameter  that  yields  the  corresponding  maximum  mean  on-axis  irra- 
diance  in  the  presence  of  jitter.  Previously,  optical  systems  modelers  have  often  had  to  resort  to 
complex  and  time  consuming  numerical  procedures  (e.g.,  Monte  Carlo  techniques)  to  account 
for  random  jitter.  The  accurate  elementary  approximations  obtained  here  allow  analytical  results 
to  be  obtained  for  the  optimum  mean  irradiance  distribution  in  the  presence  of  random  jitter, 
which,  in  many  cases,  aids  in  facilitating  parametric  estimation  and  optimization  of  over^  system 
performance. 
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2.  Far-Field  Irradiance  Distribution  in  the  Absence  of  Jitter 


In  this  report,  we  consider  an  unobscured  circular  apeiture  of  diameter  D  that  truncates  a 
(TEMoo)  Gaussian-shaped  laser  beam  of  1/e^  intensity,  radius  (Do,  and  power  P.  In  the  absence  of 
jitter,  the  irradiance  distribution  of  wavelength  X  at  propagation  distance  R  in  the  far  field  (i.e.,  for 
propagation  distances  satisfying  R  »  ic(Do^  A  )  at  angidar  coordinate  6,  can  be  expressed  as^ 

I(0)=  I(O)G(0).  (2.1) 


where  1(0),  the  on-axis  irradiance,  is  given  by 


I(0)=  lu  Ftninc* 


I  PA 
*U=  — r 

(XR)2 


P  2[l-exp{V}]^ 

*  tninc 9  * 

ii 

A  is  the  area  of  the  aperture  (=  7cD2/4),  and  G(0),  the  normalized  beam  pattern,  is  given  by 

2  J  dxx  exp|-x^  |jo  (ux/  p) 

G(0)=  - p-jr -  , 

l-exp(V) 

where  Jo  is  the  Bessel  function  of  the  first  kind  of  order  zero, 

=  Q 

“  (X/jcD)’ 


D 


and  d  =  2(Do  is  the  beam  diameter  in  the  plane  of  the  aperture.  In  Eq.  (2.2),  the  (piantity  lu  is  the 
on-axis  irradiance  obtained  for  uniform  illumination  of  the  circular  aperture  (which  is  the  abso- 


lute  maximum  on-axis  inadiance  obtained  imder  any  physical  circumstances),  and  the  quantity 
I’tninc  gives  the  effects  of  truncation  of  the  Gaussian  beam  by  the  aperture.  For  an  iqjerture  of 
fixed  diameter  D,  it  is  well  known^  that  the  maximum  on-axis,  far-field  irradiance  is  obtained  for 
d  «  0.89D.  For  this  value  of  d,  Ftnmc  0.81.  Thus,  in  the  absence  of  jitter,  the  maximum  far- 
field  intensity  is  about  81%  of  what  would  be  obtained  if  the  same  total  power  is  uniformly  dis- 
trilHited  over  the  circular  aperture;  this  condition  occurs  for  d/D  »  0.89. 

As  is  shown  in  Section  3,  the  optimum  value  of  d/D  in  the  presence  of  jitter  is  obtained  for 
(d/D)opt  <  0.89.  As  such,  we  first  derive  an  accurate  expression  for  the  beam  pattern  within  the 
central  lobe  for  d/D  <  1. 

Consider  a  beam  of  Gaussian  shape  of  the  form 

lA(0)  =  I(O)GA(e),  (2.8) 


where  1(0),  the  exact  on-axis,  far-field  irradiance,  is  given  by  Eq.  (2.2),  and 

Ga(0)=  exp[^-|(e2/0|)j, 

where 


(2.9) 


(2.10) 


and  the  parameter  b  is  determined  from  the  requirement  that  the  integrated  far-field  irradiance  of 
Ga  equds  the  exact  value  of  the  transmitted  power,  Pp,  of  the  truncated  beam.  We  have 

Pt  =  [1  -  exp(-2ji2)]p,  (2.11) 


where  p  =  D/d,  and  P  is  the  incident  power.  Thus,  b  is  obtained  as 


b  = 


(2.12) 


For  d  «  D  (i.e.,  p  »  1),  b  «  0  and  0b  =  X/Jcd,  as  expected  physically.  That  is,  for  d  «  D,  the 
expression  for  the  approximate  beam  width  reduces  to  that  of  an  infinite  Gaussian  beam.  On  the 
other  hand,  for  d  »  D,  b  »  ^^2.  That  is,  in  this  limit  the  Gaussian  width  equals  ^^2  X/nD.  Qearly, 
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for  d  >  D,  it  is  expected  physically  that  the  truncated  beam  in  the  far  held  behaves  more  like  an 
Airy  pattern  than  a  Gaussian.  However,  here  we  are  only  interested  in  beams  where  d/D  <  1.  For 
this  case,  we  have  compared,  in  Figures  1-4,  the  exact  [based  on  numerical  integration  of  Eq. 
(2.5)]  and  the  approximate  far-held  pattern  given  by  Eq.  (2.9).  Examination  of  these  hgures 
reveals  that,  within  the  hrst  central  lobe  of  the  far-held  pattern,  lA(d)  [  given  by  Eq.  (2.8)]  is  a 
very  good  approximation  to  the  actual  far-held  irradiance  pattern  of  a  truncated  Gaussian  beam. 
The  small  differences  indicated  in  the  hgures  are  well  within  the  practical  uncertainties  that  always 
exist  for  “real”  Gaussian  beams.  We  note  that  Ia(6)  gives  both  the  correct  on-axis  irradiance 
and  integrated  power  for  the  truncated  beam.  Additionally,  for  d/D  <  1,  the  integrated  far-held 
power  in  the  hrst  central  lobe  is  >  95.5%,^  and  thus  Ia  is  an  accurate  representation  of  the  beam 
in  this  regime. 

Note,  the  characteristic  1/e^  beam  half-width  angle,  Ob.  can  be  written  as 


0B  = 


2X 

JtDeff  ’ 


(2.13) 


where 


Deff  = 


i^tahh^D 


^/2d^). 


(2.14) 


Figure  1.  Far-held  irradiance  beam  pattern  of  a  truncated  Gaussian  beam  as  function  of 

normalized  angular  coordinate  u  [=  6/(X/nD)]  for  d/D  =  0.9.  The  dashed  and  solid 
curves  are  based  on  the  exact  [i.e.,  Eq.  (2.S)]  and  ^proximate  [i.e.,  Eq.  (2.9)] 
expressions  for  the  beam  pattern,  respectively. 
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Figure  2.  Far-fieldiiTadiance  beam  pattern  ofa  truncated  Gaussian  beam  as  function  of 

nonnalized  angular  cocudinate  u  [=d/(X/icD)]  for  d/Ds  0.8.  The  dashed  and  solid 
curves  are  based  on  die  exact  [i.e.,  Eq.  (2  J)]  and  approximate  [i.e.,  Eq.  (2.9)] 
expressions  for  the  beam  pattern,  respectively. 


Figure  3.  Rir-fiield  iiradiance  beam  pattern  of  a  truncated  Gaussian  beam  as  function  of 

normalized  angular  coordinate  u  [=  0/(X/nD)l  for  d/D  s  0.7.  The  dashed  and  solid 
curves  are  based  on  the  exact  [i.e.,  Eq.  (2.5)]  and  approximate  [i.e.,  Eq.  (2.9)] 
expressions  fm*  the  beam  pattern,  respectively. 
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Figure  4,  Far-fieki  irradiance  beam  pattern  of  a  truncated  Gaussian  beam  as  function  of 

normalized  angular  coordinate  u  [s  8/(X/}cD)]  for  d/D  »  0.6.  The  dashed  and  solid 
curves  are  based  on  the  exact  [i.e.,  Eq.  (23)]  and  approximate  [i.e.,  £q.  (2.9)] 
expressions  for  the  beam  pattern,  respectively. 


The  quantity  D^ff  is  the  l/e^  beam  intensity  diameter  in  the  aperture  plane  of  an  “effective” 
infuiite  Gaussian  field  that  results  in  the  distribution  Ia.  That  is,  we  replace  the  original  Gaussian 
field  of  1/e^  diameter  d  and  the  aperture  of  diameter  D  by  an  infinite  Gaussian  field  of  1/e^ 
diameter  Deff,  which,  when  propagated  to  the  far  field,  results  in  an  irradiance  distribution,  1a,  that 
is  a  very  good  s^proximation  to  what  would  be  obtained  from  the  truncated  Gaussian.  We  note 
that  Ia  does  not  correctly  describe  the  side  lobes  and  should  not  be  used  when  the  corresponding 
distribution  is  important;  it  is  a  valid  approximation  within  the  main  central  lobe  only.  In  the 
next  section,  we  use  this  infinite  Gaussian  field  to  obtain  the  corresponding  results  for  the  mean 
irradiance  distribution  in  the  presence  of  random  angular  jitter. 
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3.  Far-Field  Irradiance  Distribution  in  the  Presence  of  Jitter 


Here,  we  assume  circularly  symmetric  Gaussian  angular  jitter  beam  statistics.  Let  cp-  denote  the 
1-axis,  1-0  jitter  variance  (for  adaptive  optics  systems,  is  the  corresponding  residual  jitter  vari¬ 
ance).  Then,  assuming  a  Gaussian-shaped  beam  of  beam  width  6b  [i.e.,  Eq.  (2.9)],  it  can  be 
shown^  that  the  average  on-axis  irradiance  in  the  presence  of  jitter  is  given  by 


<I(0)>  =  I(0)Fjitter  -  luPinincPjitter. 


(3.1) 


where 


Fjitter  - 


\+-j- 

e| 


(3.2) 


is  a  reduction  factor  accounting  for  the  effects  of  jitter  on  the  mean  on-axis  far-field  intensity, 
and  angular  brackets  represent  the  ensemble  average.  In  addition,  the  corresponding 
(approximate)  mean  beam  profile  in  the  presence  of  jitter,  Ga)  .  is  given  by 


GAj(0)  =  exd 


0^ 

2 

(3.3) 


For  a  fixed  aperture  diameter  D,  examination  of  Eq.  (3.1)  reveals  that  the  dependence  of  the 
mean,  far-field,  on-axis  intensity  is  contained  in  the  term  FtmncFjitter-  Now,  Ftnmc  reaches  its 
maximum  for  d  »  0.89D.  On  the  other  hand,  Fjitter  is  a  decreasing  function  of  d.  In  particular, 
Fjitter(<l  <  0.89D)  >  Fjitter((l  >  0.89D).  Hence,  it  follows  that  the  value  of  d  that  maximizes  the 
on-axis,  far-field  intensity  in  the  presence  of  jitter  will  be  obtained  for  d  <  0.89D.  This  property 
is  shown  in  Figure  5,  where  we  have  plotted  GAj  (0)  as  function  of  d/D  for  various  values  of  nor¬ 
malized  jitter  p  [=  Oj/(X/JcD)]. 

Based  on  Eq.  (3.1),  we  can  readily  obtain  results  for  the  optimum  ratio  of  d/D  that  yields  the 
maximum  on-axis,  far-field  (or  focal-plane)  irradiance.  Differentiating  Eq.  (3.1)  with  respect  to 
“d”  and  setting  the  result  equal  to  zero  yields  a  transcendental  equation  for  the  optimum  values 
of  d/D: 


(3.4) 
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d/D 


Figure  5.  Far-fidd,  on-axis  inadiance  as  a  function  of  d/D  for  various  values  of  the  nonnalized 

jitter  standard  deviation  p  [=  aj/(X/id>),  where  aj  is  the  1-axis,  l-a  jitter  standard 
deviationl. 


where  m,  =  (D/d)opt,  and 

p  =  Oj/(X/jcD)  (3.5) 

is  the  nonnalized  standard  deviation  of  jitter. 

Althou^  Eq.  (3.2)  cannot  be  solved  analytically,  numerical  results  are  readily  obtained  and  are 
plotted  in  Figure  6  as  a  function  of  noimtdized  jitter.  Of  course,  for  p  =  0,  the  results  shown  in 
Figure  6  are  identical  to  the  results  given  in  Chap.  18  of  Ref.  1.  It  can  be  shown  that  for  p  »  1, 
(d/D)opt  l/nn(p)]^^.  An  accurate  analytical  approximation  for  (d/D)opt  for  0  ^  p  <  10,  as 
obtained  by  a  least-squares  fit  to  the  calculated  values,  is  given  by 


(d/D)opt  - 


0.892  -  0.0536p^ + 0.0089p^  for  0  ^  p  <  1 .3 


0.468-1- 


0.592  0.085 

Vp 


for  1.3  <10. 


(3.6) 


The  accuracy  of  this  approximation  is  better  than  1%  over  the  range  0  ^  p  <  10. 
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As  a  practical  matter,  it  is  expected  that  viable  optical  systems  should  be  designed  to  operate  in 
the  regime  of  p  <  2-3  or,  otherwise,  serious  degradation  of  system  performance  by  jitter  will  be 
obtained.  Over  this  range,  we  obtain  that  0.89  >  (d/D)opt  >  0.73.  For  example,  for  p  =  1  and  2, 
examination  of  Figure  6  reveals  that  (d/D)opt  ®  0.85  and  0.78,  respectively.  Figure  7  is  a  plot  of 
the  optimum  on-axis  intensity  as  a  function  of  normalized  jitter.  For  the  previous  example, 
examination  of  Figure  7  reveals  that  for  p  =  1  and  2,  the  maximum  on-axis  intensity  is  about  57 
and  30%  of  what  would  be  obtained  for  uniform  illumination  of  the  aperture,  respectively. 


Another  quantity  of  interest  is  the  encircled  trarrsmitted  power  distribution.  Denoting  by 
AP(6o)/Pt  the  lection  of  transmitted  power  contained  within  a  circle  of  radius  Oq,  we  have 


AP(eo)_ 

Pt 


00 

f  iA(e)2jted6 


0 


=  l-exp 


2efj’ 


(3.7) 


which  is  plotted  in  Figure  8  for  various  values  of  normalized  jitter,  p.  For  purposes  of  compari¬ 
son  with  the  exact  results,  which  are  based  on  numerical  integration  of  Eq.  (2.5),  we  consider,  as  a 
worst  case  of  interest  here,  the  case  d/D  =  0.89  and  obtain  that  the  encircled  transmitted  power 
contained  within  the  first  central  lobe  is  equal  to  97.0%  as  compared  with  99.2%,  as  obtained 
from  the  aj^roximation  given  by  Eq.  (3.7)  for  p  =  0.  Thus,  Eq.  (3.7)  represents  a  highly  accu¬ 
rate  representation  for  the  encircled  far-held  (or  focal  plane)  transmitted  power  distribution  for 
all  cases  of  interest  here. 


Figure  6.  The  optimum  ratio  of  the  beam  diameter  in  the  aperture  to  the 
truncation  diameter  as  a  function  of  normalized  jitter. 
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Encircled  power  distribution  Maximum  on-axis  intensity 


P 


Figure  7.  The  maximum  mean  far-fmid  on-axis  inadiance  (Le.,  the  value  few 
the  optimum  ratio  of  d/D)  as  a  function  of  normalized  jitter. 


Figure  8.  Encircled  power  distribution,  AP/Pj,  as  a  function  of  normalized  angular 

coordinate  uo  for  various  values  of  nonnalized  jitter  p.  For  each  value  of  p,  the 
beam  is  truncated  at  its  optimum  value. 
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4.  Conclusions 


In  the  presence  of  random  angular  jitter,  we  have  derived  accurate  analytical  expressions  for  both 
die  mean  far-field  irradiance  distribution  and  encircled  power  for  a  truncated  Gaussian  beam  that 
are  valid  for  d/D  <  1,  where  d  is  the  beam  diameter  in  the  aperture  plane,  and  D  is  the  truncation 
diameter.  In  particular,  the  mean  irradiance  distribution  gives  both  the  exact  on-axis  mean  irradi- 
ance,  total  integrated  transmitted  power,  and  encircled  power  distribution,  valid  to  better  than 
about  2%.  Over  the  range  of  interest,  the  difference  between  the  exact  and  approximate  results 
are  well  within  the  uncertainties  that  always  exist  in  practice.  Based  on  the  approximate  irradiance 
distribution  given  by  Eq.  (3.2),  we  then  obtained  results  for  both  the  optimum  ratio  of  d/D  that 
results  in  the  maximum  mean  on-axis  irradiance  and  the  corresponding  on-axis  irradiance.  As 
such,  the  results  obtained  here  will,  in  many  cases,  be  useful  in  facilitating  parametric  estimation 
and  optimization  of  overall  system  design. 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation's  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientiflc  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provid^  by  these  individual  Technology 
Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  in&ared  detector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optic^  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battoy  electrochemistry,  tottery  testing  and  evaluation. 

Mechanics  and  Materials  Technolt^  Centen  Evaluation  and  characterization  of 
new  materials:  metals,  alloys,  ceramics,  polymers  and  their  composites,  and  new 
forms  of  carbon;  development  and  analysis  of  thin  films  and  deposition  techniques; 
nondestructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  development  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch 
vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  and 
vulnerability  assessment;  contamination,  thermal  and  structural  control;  high 
temperature  thermomechanics,  gas  kinetics  and  radiation;  lubrication  and  surface 
phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  eff«:ts  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-field- 
of-view  rejection. 


